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ABSTRACT 

The Ap star HD 3980 appears to be a promising roAp candidate based on its funda- 
mental parameters, leading us to search for rapid pulsations with the VLT UV- Visual 
Echclle Spectrograph (UVES). A precise Hipparcos parallax and estimated temper- 
ature of 8100 K place HD3980 in the middle of the theoretical instability strip for 
rapidly oscillating Ap stars, about halfway through its main sequence evolution stage. 
The star has a strong, variable magnetic field, as is typical of the cool magnetic Ap 
stars. Dipole model parameters were determined from VLT observations using FORS1. 
From Dopplcr shift measurements for individual spectral lines of rare earth elements 
and the Ha line core, we find no pulsations above 20 — 30ms -1 . This result is corrob- 
orated by inspection of lines of several other chemical elements, as well as with cross- 
correlation for long spectral regions with the average spectrum as a template. Abun- 
dances of chemical elements were determined and show larger than solar abundances 
of rare earth elements. Further, ionisation disequilibria for the first two ionised states 
of Nd and Pr are detected. We also find that the star has a strong overabundance of 
manganese, which is typical for much hotter HgMn and other Bp stars. Line profile 
variability with the rotation period was detected for the majority of chemical species. 

Key words: Stars: oscillations - stars: variables - stars: individual (HD 3980) - stars: 
magnetic. 



1 INTRODUCTION 

After the discovery of low amplitude pulsation for the bright 
and well-stu died cool, magnetic chemic a lly peculiar A (Ap ) 
star /3CrB (jHatzes fc Mkrtichianll2004l : iKurtz et al.ll2007h . 
the question arose whether all cool Ap stars (with effec- 
tive temperatures below T e g ~ 8200 K) are rapid oscilla- 
tors. This still unanswered question could be a pivotal point 
for theoretical modelling and understanding of the driv- 
ing mechanism in rapidly oscillating (roAp) stars. Several 
searches for rapid oscillations in Ap stars have been made. 
lElkin, Kurtz, fc Mathvsl (|2008ah found that 24 known roAp 
stars that exhibit pulsations photometrically also show rapid 
radial velocity variations for the corresponding pulsation pe- 
riods. However, several Ap stars with photometric indices 
typical for know n roAp stars were tested photometrically 
for pulsations bv lMartinez fc Kurtzl (|l994h and found to be 



* Based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile, as part of programme 077.D-0150(A) 
and part of programmes 074.D-0392(A) and 076.D-0535(A) in 
the ESO Archive. 



stable to high precision; these stars are called non-oscillating 
Ap stars, or noAp stars. 

About a decade later, lElkin et al.l (|2005al ) used fast 
UVES spectroscopy to discover low amplitude radial ve- 
locity pulsation for one of these stars, the form er noAp 
star HD 116114. lLorenz. Handler, fc Kurtzl (|2005t ) reported 
a new photometric null result for the star, but noted a low 
amplitude pulsation peak in the amplitude spectrum of their 
data slightly above the noise level for the same frequency. 
Spectroscopy therefore is superior to photometry as a tool 
for detection of rapid low amplitude pulsation, and we ex- 
pect that more Ap stars previously identified as non- oscil- 
lating (noAp) may exhibit rapid radial velocity variations. 
One such promising example is HD 965, which has physical 
parameters corresp onding to tho s e of k nown roAp stars. The 
fast photometry bv lKurtz et al.l (l2003f) an d high time reso- 
lution spectroscopy by lElkin et al. (2005b), however, found 
no rapid variability. Both of these studies mentioned the 
known roAp problem of amplitude modulation with the ro- 
tation phase as a possible explanation for the null result and 
proposed to re-observe the star at a different rotation phase 
when one of the magnetic poles is at a more favourable as- 
pect. 
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As seen in, e .g., the astrometric HR-diagram by 
iHubrig et all (|2005l . their figure 2) for roAp and noAp stars, 
the apparent noAp stars occupy essentially the same regions 
as the roAp stars. However, the noAp stars app ear to be sys- 
tematically more evolved than the roAp stars l|North et al.l 
Il997l : lHan'dler fc Paunzenlfl999l : IHubrig et al.ll200d). Never- 
theless, the theoretical roAp instability strip (|Cunhall2002l ) 
also predicts rapid pulsations for the more evolved Ap 
stars (near the term inal age the main sequence). HD 116114 
l|Elkin et al.1 l2005bl ) was indeed detected in this region of 
the HR-diagram, oscillating with the predicted frequency 
0.79 mHz (the lowest frequency known for the roAp stars). 
Still the only case of a luminous roAp star, HD 116114 shows 
extremely low radial velocity pulsation amplitude and only 
for small number of chemical elements such as europium and 
lanthanum. 

Using the sam e instrument (UVES) and procedure as 
lElkin et al.1 (|2005bh . iFrevhammer et al.l (|200Sl) searched for 
rapid pulsations among a group of nine evolved Ap stars 
inside the roAp instability region in the HR-diagram. Sur- 
prisingly they did not detect any radial velocity pulsation for 
these stars, but showed that only 3 — 5 of the stars may have 
magnetic field strengths considerably in excess of 2 kG. More 
evolved stars are theoretically expected to require relatively 
strong magnetic fields to suppress local surface convection 
and facilitat e observable amplitudes of rapid pulsations (see 
ICunhall2002r ). 

The bright Ap star HD3980 (£Phe, V = 5.719 mag) 
was proposed by one of us (Hubrig) to be a good roAp can- 
didate because of its many properties in common with those 
of known roAp stars, such as a magnetic field and peculiar 
abundances of rare earth elements, as we show in Fig.[T] A 
photometric search for rapi d pulsations had a l ready been 
performed for HD 3980 by iMartinez fc Kurtj |l994|) who 
failed to detect any pulsational variability. Base d on our suc- 
cess with high time resolution spectr oscopy (e.g. lKurtz et al.l 
l2007l ; lMathvs, Kurtz, fc Elkinll2007l ). we then decided to use 
fast spectroscopy to test this promising roAp candidate for 
rapid pulsations. The following sections discuss the collec- 
tion of data, their analysis and the null result for rapid os- 
cillations. 



2 OBSERVATIONS AND DATA REDUCTION 

HD3980 was observed twice in 2006 July with UVES on 
the ESO Very Large Telescope (VLT) to search for rapid 
pulsations at two rotation phases roughly coinciding with 
the predicted times of magnetic maximum and minimum 
(see Sect. [3}, when pulsation amplitudes are expected to be 
highest. On the first observing night, 2006 July 23, 104 spec- 
tra were obtained in 2.0 hr, and on 2006 July 25, another 55 
spectra were collected in 1.0 hr. The seeing was relatively 
poor, especially on the first night when it varied from 1"5 
to 3". However, the star is bright and we used an image 
sheer (IS#3) to optimally utilise the observing conditions 
at the instrument's maximum resolution with a 0'.'3 slit. Ex- 
posure times of 40 s were then used, which together with 
readout and overhead times of ~ 28 s provided a time res- 
olution of 68 s. We used the RED (600 nm) setting which 
covers the wavelength region AA4970 — 7010 A, with a gap 
in the region AA 5963 - 6032 A caused by the space between 
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Figure 1. A small spectral region of HD3980, compared with 
those of the two roAp stars, HD 128898 (a Cir) and HD 176232 
(lOAql). The spectra of the latter two stars are shifted up and 
down, respectively, in intensity by 0.2. Selected atomic lines are 
indicated. HD 3980 shows strong lines for the rare earth elements 
Nd and Pr, a common characteristic of known roAp stars. 

the two CCD mosaic halves. The average spectral resolution 
is about R = 110 000. The camera uses two 2K x 4K CCDs 
with 15 11m pixels. Raw CCD frames were processed using 
the UVES pipeline to extract and merge the echelle orders to 
ID spectra that were normalised to the continuum. The av- 
erage signal-to-noise ratio in the continuum, estimated from 
ID spectra, is about 250. 



3 THE STELLAR PARAMETERS 

The rotation period was determined from HD 3980's d ouble 
wave light curve by iMaitzen. Weiss, fc Wool (|l980l) . who 
provided the elements (for the primary minimum): JD(Prim. 
Min. v) = 244 2314? 48 + (3?9516± d 0003) E. These authors 
also obtained 8 longitudinal magnetic field measurements 
which show large scatter when phased with the rotation pe- 
riod. For the rotation phases 0.30 — 0.32, 3 measurements 
give (B\) — —780 ± 700 G, while 4 measurements at phases 
0.78 — 0.83 give 860 ± 996 G. A single point measured near 
phase zero, with a conservative error estimated from the 
other measurements, gives 1670 ± 1000 G. These 8 measure- 
me nts are insuffici e nt to fit the magnetic field curve properly. 

IHubrig et al.l (|2006l ) obtained 3 precise measurements 
of the longitudinal magnetic field: (B\) — 1210 ± 32 G, 
395 ± 26 G and 452 ± 15 G. Two more observations were 
obtained by Hubrig, using the same method: (Bi) — — 828 ± 
18 G (JD245 4432.51) and 1804 ± 15 G (JD245 4433.54). 
These measurements are, togeth er with averaged data from 
IMaitzen. Weiss, fc Wood! (|l980l '>. shown in Fig.[2]with error 
bars and phased with the rotation period. A least squares 
sine curve fit of the five Hubrig measurements provides a 
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mean field of (Si) = 40 ± 18 G and amplitude of Ab, = 
1918 ±29 G. By intention, our two observing nights occurred 
near the times of magnetic extrema at rotation phases 0.94 
and 0.46, respectively, which are close to the positive ex- 
tremum and the negative extremum of the magnetic field 
curve; 

iHubrig. North, fc Schollerl l|2007h determined stellar pa- 
rameters: T eff = 8240 ± 310 K, logy = 4.05 ± 0.09, 
logL/L0 = 1.296 ± 0.052 and a projected rotational ve- 
locity v sin i = 15kms _1 . For the 3.95 d rotation period 
and the longitudinal field measurements, they estimated the 
magnetic field geometry parameters: Mean longitudinal field 
(Bi) = 1200 G, inclination angle between rotation axis and 
line-of-sight i = 32° and angle between rotation axis and 
magnetic dipole axis [3 = 88°. These values, however, are 
highly uncertain. 

Our fit to the magnetic measurements shown in Fig. [2] 
can also be used to constrain the inclinations of the magnetic 
and rotation axes. It is easy to show for a centred dipolar 
magnetic field that 



B\ oc Bp cos a 
and 

cos a = cos i cos (3 + sin i sin f3 cos Qt 



(1) 



(2) 



where a in the angle between the magnetic pole and the line- 
of-sight, i is the rotational inclination, (3 is the angle between 
the rotation axis and the magnetic axis, Q is the rotation 
frequency, B\ is the longitudinal magnetic field strength and 
_B P is the polar field strength. 

It is obvious from Eq. 2 that the mean magnetic field 
strength 

(B\) oc cos i cos/3 (3) 
and the amplitude of the magnetic field variations 
Ab 1 oc sin i sin (3 . (4) 
Thus we get 

tanitan/3 = Ab 1 /(B\) , (5) 

from which f3 can be constrained when i is known, us- 
ing the values for Ab 1 and (B\) from our fit to the mag- 
netic data shown in Fig. [5] We derive in Sect. 6 below 
from the Hipparcos parallax and our estimate of T e ff that 
R = 2.19jlo ig R©> which coupled to our measurement of 
vsini = 21.0 ± 3.0 km s" 1 gives a weak constraint on the 
rotational inclination of i = 49° l^. Eq. 5 can then be used 
to estimate the magnetic obliquity to be f3 = 88?6l4 g (ac- 
counting for all 1-ct errors, while it becomes undetermined 
for 2-cr errors). Several more precise magnetic measurements 
at different rotation phases are still needed to fill the gaps 
in the magnetic field curve. The important conclusion here 
is that the obliquity of the magnetic field must be near to 
90° to account for the mean magnetic field strength being 
close to zero. 

The stellar parameters T ff and log o were estimated 
from Stromgren photomet ry by iMartinea (|l993) using the 
iMoon fc Dworetskvl (|l985h calibrations. Then a small grid 
of synthetic Ha line profiles was compared to the ob- 
served spectra. The model spectra were produced using Ku- 
rucz model atmo spheres (|Kuruczl I1979T ), models from the 
NEMO database (|Heiter et all 12002^ and calculations with 



the SYNTH program by IPiskunovl (|1992T ). It is difficult to 
locate accurately the continuum in the broad Ha region. 
The resulting normalised profiles are slightly different for 
the two observing sets. This may partly be due to our con- 
tinuum placement, but the Ha profile may itself be variable. 
An exam ple of this was found in t he ho tter peculiar star 
36 Lyn by iTakada- Hidai fc Aikmanl (|l989l ). Variation in the 
Ha profile is consistent with small variations in the effec- 
tive temperature of the star as it is rotates. For the average 
spectrum of the first observing night, we obtained a good 
fit for T c ff = 8000 K and logp = 4.0, whereas for the second 
night T c ff = 8200 K was required. 

We performed an abundance analysis of the first set 
of spectra of the star, for which T c ff = 8000 ± 200 K and 
logp = 4.0 ± 0.2. A Kurucz model atmosphere with T c ff = 
8000 K and logo = 4.0 and with enhanced (above solar) 
metallicity ([M/H]= +0.5) was used for calculating model 
spectra with SYNTH. Spectral line lists were taken from the 
Vienna Atomic Line Database (VALD, Kupka et al.l 1999T) 
and the DREAM database (|Biemont. Palmeri. fc Quinetl 
Il999ft . We determined vsini = 21.0 ± 3.0kms _1 from sev- 
eral symmetric lines. This val ue is slightly higher than 
the v sin % = 15kms -1 found by IHubrig. North, fc Schollerl 
(|2007l ). That could be explained by, e.g., a variable mag- 
netic field strength which contributes to the broadening of 
lines with large Lande factors. A more probable explanation 
is line profile variations due to a spotted abundance dis- 
tribution. Some spectral lines are indeed particularly nar- 
row and give a lower value of vsini. We consider this to 
be an effect of these elements being concentrated in spots 
on the stellar surface such that their spectral lines do not 
reflect the full rotational broadening. An example is Li I 
for which spot ted distributions have been demo n strate d for 
roAp stars by Faraggia na. Gerbaldi. fc Delmasl (|l996T l and 
IPolosukhina et al.l ( 2000h . For the strong, isolated and sym- 
metric but non- resolved doublet Li I 6707.76 and 6707.91 A 
we obtained vsini = 16.5kms _1 . 



4 SEARCH FOR RAPID RADIAL VELOCITY 
VARIATIONS 

We searched for pulsations as periodic Doppler shifts in two 
ways: first by using cross-correlations of long stretches of 
spectral regions, then by measuring centre-of-gravity shifts 
of individual lines in the spectra. For frequency analyses, 
we used MIDAS's TSA (Time Series An alysis) conte xt, a 
discrete Fourier transform programm e bv lKurtzll|l985h and 
the PERIOD04 (|Lenz fc Bregerll2005r i programme. 



4.1 Cross-correlation radial velocity analyses 

The cross-correlation method, using large spectral regions, 
is often useful for detecting pulsation in roAp star candi- 
dates and for finding additional frequencies in known roAp 
stars (see, e.g.. iMatthews et al.lll988l : iBalona fc Zimall2002l : 



lHatzes fc Mkrtichianl l2004h . The cross-correlation ampli- 
tudes from correlation of long spectral regions are, though, 
not directly comparable to those derived from line pro- 
file measurements. This is mainly due to different pulsa- 
tion amplitudes and phases of different ions in the strati- 
fied roAp atmospheres, where low amplitude elements such 



4 Elkin et al. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Rotation phase 



Figure 2. The longitudinal field observations available for 
HP 3980. T he five filled circles are measurements from 
iHubrig et alj d2006h and two recent points observed on 
JD245 4432.51 and JD245 4433.54 (Hubrig, unpublished). The 
full line shows the least squares fit of a sine function to these 
5 points, using phases based o n the light curve elements by 
iMaitzen. Weiss, fe Woodl <198Ch. For compari s on, 3 measurc- 
ments (triangles) by llVlaitzen. Weiss, fe Woodl (ll980h are also 
shown with standard deviation error bars and a weighted fit 
(dashed line) to all 8 points. Good agreement with the recent 
and old data is seen, in spite of the long time baseline. 



as Fe dilute the 'integrated' Dop pler shifts. An example 
is the r oAp star 10 Aql for which lElkin. Kurtz, fe Mathvsl 
(|2008bh found a pulsation amplitude of 6 ± 1ms 1 with 
cross-correlation for the spectral region 5150-5450 A, while 
individual lines of rare earth elements showed amplitudes in 
excess of 500 ms -1 . 10 Aql represents a case of a roAp star 
for which the lines of rare earth ions are weak. Another 
example is HD 154708 which has a very strong magnetic 
field. From cross-correlation mea surements of HD 154708 in 
the spectral range 5150 - 5800 A iFrevhammer et~aT1 (|2008u 
detected pu l sation with an amplitude of 10 ms -1 , while 
iKurtz et all (|2006al ) obtained amplitudes around 60 m s 
for some individual rare earth lines. 

Cross-correlations were performed with our HD 3980 
spectra, using the average spectrum as template. For the 
line-rich spectral range 5000 - 5960 A, we searched the re- 
gions 5150 - 5400 A and 5400 - 5700 A, but in both cases 
no significant signal was detected above a level of 10 ms -1 
with a = 3ms -1 (see Fig. [3] top panel) in the frequency do- 
main of known roAp stars. Only low frequency peaks below 
0.4 mHz are seen. However, we disregard these periodicities 
as they are also seen in the radial velocity measurements 
of telluric lines, hence are instrumental in origin. A spec- 
tral region from 6863 — 6938 A with abundant telluric lines 
was used to check the instrumental stability and identify 
non-stellar periodicities. The results from telluric lines only 
show low frequency noise (due to instrumental drifts or me- 
teorological changes). The telluric region otherwise shows 
stability at the level below 10ms -1 with a = 3ms -1 . The 
spectral region longwards of 6000 A has a lower line density 
and a higher scatter, but confirms these null results. 



4.2 Line profile radial velocity analyses 

In roAp stars, lines of rare earth elements typically show 
the largest Doppler shift pulsation amplitudes. Amplitudes 
vary for different elements and range from a few dozen me- 
tres per second up to a few kilometres per seco nd for var- 
ious roAp stars (|Elkin. Kurtz, fe Mathvsl 12008 j ) . Also the 
narrow line co re of the Ha profil e shows rapid pulsations 
in roA p stars jKurtz et all l2006al ; lElkin. Kurtz, fe Mathvsl 
l2008al ). We therefore searched for pulsations in HD3980 us- 
ing the centre of gravity method for similar lines and show 
amplitude spectra of selected lines in Fig. [3] Although we 
concentrate on analyses of lines of the rare earths, other 
chemical elements were also tested. All radial velocity curves 
subjected to period searches were de-trended beforehand 
with linear least square fitting, which also eliminated the 
barycentric velocity correction which is approximately lin- 
ear during the short duration of our two runs. The Ha core 
was stable above the 20ms -1 with a noise level of 7ms -1 . 

Of the rare earth elements, we examined lines of Cell, 
Prill, Ndni, Ndn, Eun and Gdn, but no pulsations were 
detected from any of these lines above typical upper limits 
of 30 ms -1 (with noise levels varying from 10 to 30 ms -1 
for the majority of good lines). Combination of 3 — 4 lines 
reduced the noise level of individual lines, but also did not 
show any reliable pulsation signal. Due to the rotational 
broadening and many asymmetric line profiles, the number 
of suitable strong, symmetric lines was rather limited. Of the 
non-rare-earth element lines analysed, including the strong 
sodium D and Mgl lines and several Ca, Sc, Ti, Cr, Fe and 
Ba lines, no pulsations were detected above 15 — 30 m s -1 . In 
general, some lines show intriguing peaks in the amplitude 
spectra for the frequency range typical for roAp stars. But 
as these peaks could not be confirmed by other spectral lines 
of the same element or by other rare earth lines, they were 
rejected as indications of pulsation in HD 3980. Examples of 
such peaks are seen in Fig.[3]for the asymmetric and blended 
line of Pr ill 6090 A which shows two peaks with frequencies 
0.546 mHz and 2.752 mHz. Combining radial velocity mea- 
surements of three other Pr ill lines (seen in the panel below 
that for Prill 6090 A) shows no significant peaks for these 
two frequencies nor for any others. 

Among the Ap stars, HD 3980 has strong lines of lithium 
(e.g. equivalent width of 101mA for Li I 6708 A). We there- 
fore also searched the Li I 6708 A doublet for pulsation. This 
line shows strong variability with rotation phase, but is sta- 
ble to rapid pulsations above 30 ms -1 with a noise level of 
11ms -1 . 



4.3 Linear trends in line profile radial velocity 
shifts 

Line shapes of Nd, Eu and Li are observed to change strongly 
with rotation phase in HD 3980, such as the example in Fig. [3] 
with Ndni 6145 A. This is a clear indication of an inhomo- 
geneous surface distributions of elements that appear to be 
located in patches. The line analysis is complicated by the 
resulting line asymmetries, as well as by many lines that are 
blends because of HD 3980's relatively high rotation velocity 
and strong lines of peculiar elements. 

A large number of spectral lines show quasi-linear 
drifts in radial velocities which appear as low frequency 
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Figure 3. Amplitude spectra for spectral lines of HD 3980. there 
are no significant peaks. 



peaks in their amplitude spectra. Different lines may show 
different drifts over the interval of 1 — 2hr. We have 
noticed the sa me effect for several oth er Ap and roAp 
stars (see, e.g.. iFrevhammer et~aD 120081 ). For HD99563, 
lElkin. Kurtz, fc Mathvsl (|2005ch detected and discussed 
even l arger trends for particu lar lines. For the roAp star 
(5 CrB. lKurtz fc Leond (|2006bf ) similarly detected non- linear 
low frequency trends which they de-trended from the ra- 
dial velocity series. For the moderately fast rotating Ap 
stars, the trends are partly caused by spots seen at differ- 
ent aspects. For example, HD 99563 has vsini = 30kms _1 
and HD3980 has vsini = 21kms~ ; both these stars show 
stronger trends than those seen for the slower rotator f3 CrB 
(v sin i = 3.5 km s _1 ). We observe that these trends also vary 
from element to element, depending on the particular sur- 
face distribution. 

For HD 3980 we de-trended linear drifts prior to the ra- 
dial velocity analyses. However, those trends clearly contain 
physical information about the star's surface distribution. 
The Li I 6708 A doublet shows one of the strongest trends 
with a linear radial velocity change of about 0.77 kms -1 h . 
We believe that this is the result of a concentrated strong 
spot of Li seen at varying aspects with rotation. Consider- 
able line profile changes with the rotation phase are clearly 
visible in this line during our 2hr run supplemented with 
existing spectra in the ESO Archive at several other rota- 
tion phases. The available data only allow a crude tomo- 



graphic study, but do indicate different surface distributions 
of various elements. A full Doppler imaging study (such as 
iLehmann et al.ll2007l ) of this star over its rotation period is 
needed. 



4.4 The second observing run 

The second observing run, obtained about half a rotation 
period (0.49 phase difference) after the first run, collected 
55 spectra in 1 hr. Though only half the length of the first 
run, this set is sufficient for searching for frequencies in the 
known range for roAp stars. As for the previous run, this 
set also shows linear radial velocity drifts of comparable, or 
slightly higher amplitudes. Again, we see drifts for the rare 
earth element lines with increasing radial velocity, while lines 
of iron peak elements and light elements show no drifts at 
all, consistent with spots for the rare earth elements and a 
more uniform distribution for the lighter elements. 

After de-trending the radial velocity series for individ- 
ual lines, the frequency analysis again showed that there is 
no pulsation above amplitudes of 30 ms" 1 for noise levels of 
a — 10ms -1 . Cross-correlations of this set of spectra with 
their mean spectrum as template excluded periodic variabil- 
ity above 10ms _1 (a — 3ms" 1 ). 

Line profile shapes in the average spectrum are very 
similar to those of the first run, although they cover two op- 
posite sides of the spotted stellar surface at the moments of 
maximum and minimum magnetic field strengths. Spectral 
lines such as Lin and Eull that also have rather compara- 
ble line profile shapes, differ strongly from the same lines 
in spectra from the ESO archive at rotation phases 0.2 and 
0.8. For example at rotational phase 0.77, the lines of Ndm 
(see, e.g., Ndlll 6145 A in Fig. [4} have double profiles with a 
sharper, stronger and more blue component. For this phase, 
the data suggest that Nd is concentrated in two different 
spots, one at each magnetic pole. In both our observing runs 
this line seems more symmetric and wide, which may corre- 
spond to an extended surface region with overabundance of 
neodymium. 



5 CHEMICAL ABUNDANCES 

The presence of high overabundances of rare earth elements 
is one of the characteristics of known roAp stars. Another 
important property is that their abundances of the first two 
ionised states of neodymium and praseodymium show more 
than 1 dex difference, with the doubly ionised ions, which 
form higher in the atmosphere, being the most abundant 
(|Rvabchikova et al.|[2004l ; lKurtz et al.ll2007i ) . This ionisation 
disequilibrium anomaly may be explained mostly by con- 
centration of rare earth elements in high atmospheric lay- 
ers (stratification) and p artly by non-LTE (NLTE) effects 
l|Mashonkina et all [2005). NLTE effects may enhance the 
ionisation of Nd II or Pr II and accordingly weaken their ab- 
sorption lines, while strengthening those of the second ioni- 
sation state. The disequilibrium may not be limited to roAp 
stars, which are relatively cool Ap stars, but may also exist 
for hotter Ap stars. 

To compare HD 3980 spectrally with the roAp stars, we 
measured its chemical abundances and tested for ionisation 
disequilibria for Nd and Pr. Abundances were determined 
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Table 1 . Chemical abundances for HD 3980 for sc- 
lected elements, and their corres ponding solar abundances 
iAsplund, Grevesse. &: Sauvall 120051 ). The errors quoted are 
internal standard deviations for the set of lines measured. 
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Figure 4. Change of the Ndm 6145.07 A line profile with rota- 
tional phase. Our search for pulsation was done for phases corre- 
sponding to the upper and lower spectra. 



by fitting synthetic spectra to the observed average spec- 
tra. Model spectra were calculated as described in Sect. [3] 
The resulting abundances for the 2hr run average spec- 
trum are presented in TablefT] Sol ar abundances are from 
lAsplund. Grevesse. fc Sauva 3 120051 ). Only lines of Mg I sug- 
gest slightly less than solar abundances. The Fe abundance 
is near solar, or slightly greater. Notably, we detected large 
overabundances of manganese. This element typically shows 
overabundances for the hotter group of peculiar HgMn stars, 
while for cooler roA p stars, Mn mostly show s abundances 
slightly above solar (|Rvabchikova et al.ll2004l ). 



The rare earth elements and lithium (Fig. [5} show abun- 
dances that are much greater than solar. We find indica- 
tions of ionisation disequilibria for both neodymium and 
praseodymium: A[Nd]ni-n = 1-4 ± 0.3 and A[Pr]ni-n = 
1.2 ± 0.2 (where the errors are calculated from the error in 
the mean abundance for each ion's set of lines), consistent 
with th e spectroscopic signature for most roAp stars discov- 
ered bv lRvabchikova et all |2004 ). Thus, we find HD 3980 is 
spectrally comparable to known roAp stars. Three lines of 
terbium also suggest ionisation disequilibrium for this ele- 
ment. 



6 DISCUSSION 

The cool Ap star HD 3980 shows no rapid oscillations in 
radial velocity above a few tens of ms" 1 at two rotational 
phases separated by half a rotation period. We have shown 
that at least one of these coincides with a rotation phase 
near a magnetic extremum when pulsation amplitude for an 
oblique pulsator is expected to be highest. The methods ap- 
plied, and the lengths of the data sets acquired, would be suf- 
ficient for detecting pulsations in all known roAp stars. We 
therefore now reconsider whether HD 3980 is a good roAp 
candidate, when comparing its characteristics to those of 
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0.10 


Mgl 


3 


7.23 


± 
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known roAp stars. If indeed it is a good candidate, it either 
does not pulsate (thus is a noAp star), or we failed to detect 
the pulsations, or, alternatively, the roAp class character- 
istics are still too poorly established to physically discern 
noAp stars from roAp stars. 

The blue edge of the roAp instability strip where 
it crosses the mai n sequence is not firmly established. 
iHubrig et all (|2000h showed that it is ar ound 85 K, w hile 
the theoretical roAp instabi lity strip by Cunhal (l2002h ex- 
tends up to around 9500 K. iRvabchikova et al] l)2004l ) sug- 
gest a transition region (noAp/roAp) around 8100 K where 
cooler Ap stars may be rapid oscill ators. HD 3980 has a pre- 
cise parallax (n = 14.91 ± 0.35, Ivan Leeuwenl l2007h that 
for A v = 0.054 (NASA/IPAC IRSA maps) gives a lumi- 
nosity of logL/L0 = 1.27 ± 0.03 (including all uncertain- 
ties). We found an effective temperature of 8100 ± 200 K 
(estimated error) which places the star well inside the the- 
oretical roAp instability strip. The luminosity was derived 
using a bolomet ric corr e ction , BC = 0.024 ± 0.009, from 
the relations by iFlowerl (|l996T l for the range of T a g within 
its error. From this luminosity and T e g we derive a radius 
of R = 2.19t 019 
uity to be /3 

ously difficult to determine for A stars with peculiar abun- 
dances. If one similarly uses the bolometric calibration by 
lLandstreet et all (|2007T l. accounting for a 0.1 mag intrinsic 
error in the calibration, the corresponding luminosity and ra- 



2.19lo';[g R0. Eq. 5 then gives the magnetic obliq- 
?7jl4Q. Bolometric corrections are notori- 
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Fig ure 5. Relative abun dances (log e — log eq ) for HD 3980 (cir cles with error bars) co mpared to those of the roAp star /3 CrB (diamonds) 
from lKurtz et al ] l|2007t ). and the Am star 32Aqr (asterisks) bv lAdelman et al . (1997) showing the strong overabundances typical of the 
cool Ap stars. 
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Figure 6. A theoretical HR-Diagram for a sample of roAp 
stars. Luminosi ties were calculated us ing Hipparcos trigonomet- 
ric parallaxes l lPerrvman et al l |l997|) except fo r HD 3980, for 
which the revised parallax of Ivan Leeuwenl 0007) was used. The 
iMoon fe Dworetskvl jl985l) calibration of Stromgren photometric 
indices was employed for estimating the effective temperature. 
The position of HD 3980 is shown by the filled diamond which is 
the size of the error bar in logL/LQ. Solid lines ar e evolutionary 
tracks for stars with 1.5, 1.7 and 2.0 Mq taken from lSchaller et al] 
lll992l ). 



2.14 



+0.30 



-0.26 R Q 



dius become: log L/Lq = 1.25 ± 0.07 and R : 
resulting in /3 = 88?6t°8 9 4- 

The roAp stars have strong magnetic fields that range 
from several hund red Gauss to 24.5 kG in HD 154708 
(jHubrig et alj|2005l ). HD 3980 has a longitudinal field, which 
is within this range. The evolutionary tracks in Fig. [6] as well 



as those in figure 1 of lCunhal (|2002l ), indicate that HD3980 
in the midst of its life on the mai n seque n ce, ha lfway evolved 
from the ZAMS to the TAMS. ICunhal (|2002r i argued that 
more evolved roAp stars require stronger magnetic fields 
to suppress the upper envelope convection and enable the 
rapid oscillations to reach observable amplitudes. She fur- 
ther speculated that the magnetic field intensities needed 
for suppressing the convection in roAp stars more often are 
found in roAp stars with magnetically resolved lines. That 
typically requires a magnetic field modulus of ~3 kG and a 
slowly rotating Ap star (vsini = 1 — 3kms _1 ). Based on the 
longitudinal field measurements of HD3980, it could have 
a magnetic field modulus more then 3kG. The rotational 
broadening dominates any magnetic splitting of lines at this 
field intensity and direct measurement is not possible. 



iFrevhammer et all (|2008h searched for pulsations in 9 
evolved cool Ap stars located inside the theoretical instabil- 
ity strip, the majority having estimated temperatures below 
8100 K. With similar precision and upper limits on radial 
velocity amplitudes as in this study, these authors found 9 
null results. Out of 7 stars, only 3 had magnetic fields sig- 
nificantly stronger than 2kG, which possibly explains most 
of their null results for such evolved stars. The magnetic 
field of HD 3980 is strong enough to suppress local convec- 
tion and enable pulsational driving of rapid pulsations to 
reach observable amplitudes. The star does not appear to 
be near the terminal end of its main sequen ce lifetime, in 
which case this explanation may be less likely. iHubrig et al.l 
(|2000l ) pointed out the apparent deficiency of close bina- 
ries among roAp stars (although a 'handful' of exceptions 
are known). In this context, it is interesting to note that 
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a p r 

binary (see, e.g.. IPerrvman et al.lll997h . 

Our abundance analysis found that HD 3980 has strong 
enhancement of rare earth elements, ionisation disequilib- 
ria for Nd and Pr and inhomogeneous surface distributions 
of these elements. In these respects, HD 3980 strongly re- 
sembles known roAp stars. Only a relatively high temper- 
ature, compared to the typical roAp stars, may explain 
why no pulsations are found. We cannot exclude extremely 
low amplitude pulsations of about 5 — 10 ms -1 . Low am- 
plitude oscillations in ro Ap stars are known from, e.g. , 
(JCrB (|Kurtz et al]|2007tl. HP 154 708 |Kurtz et all 1200631 ) 
and HD 116114 (|Elkin et alj|2005ah . For such low pulsation 
amplitudes, ms _1 or even su b-ms -1 high prec i sion m ea- 
surements are required, such as lKochukhov et al.l (|2007l ) ob- 
tained with HARPS for HD 75445. This in turn requires the 
stars to be slow rotators, although cross-correlation for long 
spectral regions partly compensates for that. In the case of 
HD 3980, even the high precision of cross-correlation analy- 
sis did not detect low amplitude pulsations. 

It is important to find roAp stars with rotation periods 
of a few days, such as for HD 3980, as they can be subjected 
to 3D pulsation and a bundance studies, such as HR3831 
ijKochukhov et aI1l2007l ). But if the pulsation amplitudes are 
very low, such a study will be limited to very few individual 
lines or elements and will not be very informative. 

We have found that HD3980 is a bright, nearby, Li- 
spotted Ap star with a short rotation period which makes 
it an excellent observing target for Doppler imaging. We 
demonstrated that the surface distribution is spotted for 
several elements, such as Nd, Eu and Li. We have shown 
that the star has similar characteristics to the known roAp 
stars, but that it is stable to pulsations with an upper limit 
to the radial velocity amplitudes of individual spectral lines 
of 30ms _1 , less than the amplitudes of the known roAp 
stars. 
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